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Abstract 

We study the effect of migration between coupled populations, or patches, on the 
stability properties of multistrain disease dynamics. The epidemic model used in 
this work displays a Hopf bifurcation to oscillations in a single well mixed popu- 
lation. It is shown numerically that migration between two non-identical patches 
stabilizes the endemic steady state, delaying the onset of large amplitude outbreaks 
and reducing the total number of infections. This result is motivated by analyz- 
ing generic Hopf bifurcations with different frequencies and with diffusive coupling 
between them. Stabilization of the steady state is again seen, indicating that our 
observation in the full multistrain model is based on qualitative characteristics of 
the dynamics rather than on details of the disease model. 
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I. INTRODUCTION 



In this work we study the stability of a multistrain disease model in two coupled popu- 
l ations. Multistra i n dise ases ar e diseases with m ultiple coexist ing strains, such as infl uenza 



rtAndreasen et al. 



19971 ) , HIV (IHu et all ll996T ) , and dengue ([Ferguson et al. 



1999b ). We 



consider two kinds of strain interactions: cross immunity and antibody-dependent enhance- 
ment. When the disease infects an individual, his or her immune system creates serotype- 
specific antibodies, which will protect the individual against that se rotype. Howeve r, anti- 
bodies may also give some cross-protection to the other serotypes (jHalsteadl . 120071 ). This 
reduced susceptibility to the other serotypes may be temporary. When the temporary 
cross immunity wanes, heterologous secondary infections are possible. Low level antibod- 
ies developed from primary infections may fo rm complexes with t he virus so that more 
cells are affected, and viral load is increased (IVaughn et all . |2000| ). This effect is called 
antibody- dependent enhancem ent (ADE), an d it has been observed in vitro in diseases such 
as ebola ( Takada et al. . 2003 ) and dengue (jHalstead and O'Rourke . 1977 ). Throughout 
this work we make the hypothesis that AD E increases the infectiousness of secondary in - 



fective cases due to the higher viral load (ICummings et al. 



2005 



Schwartz et al. 



20051 ) . 



Alternative views o f ADE as an increase in mortality associated with secondary infectives 



can be considered (IKawaguchi et al 



20031 ) . We focus in this paper on the multistrain 



disease dengue, which is b elieved to exhib it both temporary cross immunity and antibody- 
dependent enhancement (jHalsteadl . 120071 ). Dengue is a subtropical mosquito-borne dis- 
ease that exhibits up to four serotypes. It is widespread in tropical regions of southeast 
Asia, Africa, and the Americas, infecting an estimated 50 to 100 million people every 
year (IWorld Health Organization Website! . 120061 ). Primary infections may be asymptomatic, 
while secondary infections are more severe, with about 5% of secondary infections leading 
to dengue hemorrhagic f ever (DHF) or dengue shock syndrome (DS S), the potentially fa- 
tal forms of the disease (IWorld Health Organization Website! . 120061 ) . An effective vaccine 
against dengue is very difficult to achieve. Because of ADE, infection with an unvaccinated 



strain following a single str ain vaccination may 
with secondary infections (jHalstead and Deed . 



ead t o the more severe symptoms associated 



20021 ). Therefore, an effectiv e vaccine must 



prote ct against all four serotypes simultaneously. A recent theoretical study (jBillings et al 



20081 ) has shown that eradication of dengue using only single-strain vaccines is unlikely. Be- 
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cause a tetravalent vaccine is not immediately forthcoming, deepening our understanding of 
the dynamics of dengue in more realistic models is of great importance. 

The dynamics of dengue in a single, well mi xed population has b een studied in seyeral re - 



cent publications, including but not lim i ted to 



Ferguson et al.l (jl999al ); 



Bianco et al 



Schwartz et al. 



(2009) 



Cummings et al 



(120051 ) ; IWearing and Rohanil (120061 ) . ADE and cross 



(J2005J); 



immunity have been shown to play a fundamental role in the spreadi ng of dengue , caus- 



ing instability, d e synchronization of serotype s, and chaotic outbreaks (IBianco et al 



2009 



Cummings et all 



2005 



Schwartz et al. 



20051 ). However, real populations may be spatially 



heterogeneous. To gain further insight into the wide spectrum of possible dengue dynam- 
ics, we relax the assumption of a well mixed population. The division of a population into 
spatially distinct patches simulates the potentially heterogeneous environment in which the 



disease spreads. Spatial heterogeneity has b een invoked in the past to account for t 



sistent character of some infectious diseases (IHagenaars et al. 



2004 



and to explain in-phase and out-of-phase dynamics of diseases (IHe and Stone 



r e per - 



Llovd and 



Mav 



200.T L 



1996|) 



A 



multipatch model with a ge structure has been used to explain bi- and tri-annual oscillations 
in the spread of measles (IBolker and Grenfelll . Il995l ) . 

Including human mobility patterns has a significant influence on the spreading of in- 
fectious disease. We will assume here that the coupling between p atches is via migra - 



tion, movement of individuals 



Liebovitch and Schwartzl (120041 ); 



rom one p a tch in t o another (as in 



Ruan et al. 



Grais et al 



(120031 ): 

(120061 ): ISattenspiel and Dietz! (Il995h ). An al- 



ternative coupling strategy is mass action coupling, in which susceptibles in one patch 



are assumed to interact wit 



terms (IBolker and Grenfell 



1 infectives in anothe r patch, introducing 



1995 



Lloyd and Mavl. 



1996 



Rohani et al. 



cou pling has also been modeled (see iKeeling and Rohanil (120021 ) and references therein, 



nonlin ear coupling 



19991). Stochastic 



and 



Colizza and Vespignanil (120081 )). 



The purpose of this work is to analyze the stability of a model for multistrain diseases 
with interacting strains, using dengue as an example, in a system of two coupled patches. 
Since chaotic outbreaks are likely to produce a higher number of infected individuals, under- 
standing the stability properties may play an important role in public health. The case of 
non-identical parameters in the two patches will be of particular interest, as this serves as a 
model for spatial heterogeneity. The paper is divided as follows. In Section [Til we introduce 
the epidemic model. Section [TTTI summarizes the bifurcation structure for a single patch. In 
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Section [TV] we present numerical results for bifurcations in two coupled patches. Section IVl 
motivates these results via analysis of a simple, lower dimensional model, and Section IVII 
concludes. 



II. MODEL 

We use a compartmental model for multistrain disease spread with cross immunity and 
antibody- dependent enhancem ent, which was previously studied in a single well-mixed pop- 



ulation by 



Bianco et al. 



(120091 ). In this model, individuals may develop a primary infection 
with any of the serotypes. Immediately after recovering, the individual experiences a period 
of temporary partial cross immunity to all other serotypes. When the cross immunity wears 
off, immunity to the primary infecting strain is retained, but the individual may develop 
a secondary infection with a different serotype. Infectiousness of the secondary infectives 
is increased due to antibody-dependent enhancement. After the secondary infection, com- 
plete immunity to all serotypes is assumed. A flow diagram for the single patch model with 
two serotypes is shown in Figure [T] for simplicity, but we present results here for all four 
serotypes. We extend the model to two spatially distinct patches, which are coupled by 
linear migration terms. For two patches (indexed by q) and n strains (n arbitrary), the 
model is as follows: 

~T~ = fJ'q ~ fiqSq ^ ] I Xqj + (f) ^ ] X q,ji I ~~ ^q s q 
i=l V d +i J 

-VSq + VSqi (1) 

dx r . 



<>q,% 



q,i 



II PqSq I Xq^i + (j) ^ ^ x q,ji J & x q,i flqi 

-VXq^i + VXqi^i (2) 



-OCqj - HgCq^i - VCgj + VCqi^i (3) 

= ® c q,i ~ Pq r q,i ^2 I x q,j + '^2 x Q> k i I ~~ ^q r q,i 

-VTq^i + VTq^i (4) 



-- [3 q r q>i j Xgj + x. 



dt 



l^q%q,ij ^^q,ij ^Xqi (5) 

where the variables are s g , the fraction of susceptibles in patch q; x q ^, the fraction of primary 
infectives with strain i in patch q; c q ^, the fraction of individuals in patch q with partial 
cross immunity to strain i; r q ^, the fraction of individuals in patch q that are recovered from 
a primary infection with strain i and no longer have cross immunity to the other strains; and 
Xq^j, the fraction of individuals in patch q that were previously infected with strain i and 
currently infected with strain j. The parameters are the number of strains n, the contact 
rate in patch q (3 q , the recovery rate a, the ADE factor 0, the strength of cross immunity e, 
the rate 9 for cross immunity to wear off, the birth and mortality rate in patch q n q , and 
the migration rate between patches v. A list of parameters appears in Table [H 

For simplicity, the birth and death rates in a patch are set equal to each other so that the 
total population of each patch is constant. The model of Eqs. [T]l5] allows for one reinfection. 
The parameter e determines how susceptible the cross immune compartments q are to a 
secondary infection, where e = means no cross immunity (the infection rate is identical 
for compartments q and r^) and e = 1 confers complete cross immunity (cross immunes are 
immune to a secondary infection for an average time 9~ l ). We allow e to take any value 
between and 1. The ADE factor (f> is the enhancement in the infectiousness of secondary 
infectives. 0=1 means that secondary infectives are as infectious as primary infectives, 
while = 2 m eans that secondary infe c tives are twice as infectious as primary, and so forth. 



In contrast to lAguiar and Stollenwerkl (120071 ) , we will not consider values of smaller than 
1. The migration rate from patch q to q', and from patch q' to q, is v. All individuals migrate 
with equal probability, independent of their infection status. The migration rate is assumed 
to be slow compared to the infection spread. For convenience, we put it on the same order 
as the birth/death rate. We assume that the social parameters, which are the contact rate 
(3 q and birth/death rate may vary between patches, while the epidemic parameters are 
the same in all regions. Because the social parameters depend on human factors (and in 
the case of the contact rate also include mosquito levels, which are weather-dependent), 
these parameters are the most likely to be different in adjacent regions. We use parameter 



TABLE I Parameters used in the model 

Parameter Value Reference 



/x, birth and death rate, years 
f3, transmission coefficient, years 
a, recovery rate, years -1 
9, rate to leave the cross 
immune compartment, years -1 
cfi, ADE factor 

e, strength of cross immunity 
v, migration rate, years -1 
n, number of strains 



~ 0.02 Ferguson et al. (1999) 

~ 200 Ferguson et al. (1999) 

50 Rieau-Perez et al. (1998 ) 
2 Wearing and Rohani (2006) 

> 1 Schwartz et al. (2005) 

0-1 Bianco et al. (2009) 
) - 0.05 
4 



values compatible with dengue fever, which are summarized in Table [H Our contact rate 
(3 corresponds to a reproductive rate of infection -Rn of 3.2 — 4.8, w hich is consistent with 



previous estimates ([Ferguson et al 



1999bl : iNagao and Koelle 



20081 ) 




2003 ). Note the 



FIG. 1 Flow diagram for single patch model with 2 serotypes (Bi anco et al. 
reduction of susceptibility to a secondary infection through the cross immunity factor (1 — e) and 
the enhancement of secondary infectiousness due to the ADE factor <p. Mortality terms for each 
compartment are not included in the diagram for ease of reading. 



III. SINGLE PATCH BIFURCATION STRUCTURE 



A similar model to the one of Eqs. [T}{5] has recen tly been usee 
of dengue fever in a single, well mixed, population (jBianco et al. 



;o an alyze the dynamics 



2009|). The ADE and 



cross immunity strength e were varied as bifurcation parameters. In the absence of cross im- 



munity (e 



0), APE alone generates instability, desync hronization, and ultimately chaotic 



outbreaks (jCummings et al. 



2005 



Schwartz et al 



20051 ). A Hopf bifurcation is observed for 



a critical value of the ADE factor 0, above which oscillatory solutions are obtained. Weak 
cross immunity stabilizes the syste m, while strong cros s immunity triggers instability and 



chaos even in the absence of ADE (IBianco et al. 



20091 ). In the latter case, destabilization 



occurs via a Hopf bifurcation for a critical value of the cross immunity strength e. At the 

complex pai r s of e igenvalues of the Jacobian simultaneously be- 



Bianco et al. 



(120091 ) discusses the full two parameter bifurcation 



bifurcation, three identica 
come unstable. Although 
structure (in e and 0), we will consider the cross immunity and ADE effects separately in 
the present work. 

Figured shows the bifurcation behavior of the single patch model with no cross immunity 
as the co ntact rate is varied . The bifurcation structure was computed using a continuation 



routine (IDoedel et al. 



19971 ). The ADE value at which the bifurcation occurs increases 



slightly as increases, and the dependence is approximately linear. The period of the 
periodic orbit, shown in Figure Wp for a fixed ADE value, also varie s with 0. (Note tha t 
multistrain models with ADE can display subcritical Hopf bifurcations (IBillings et all 120071 ). 
so the periodic orbit shown in Figure [2]d exists throughout the range of values shown.) 
Similarly, varying the contact rate in the absence of ADE affects the location of the Hopf 
bifurcation in cross immunity e and its characteristic frequency of oscillation. Likewise, 
varying the other social parameter, the birth rate /x, affects the location and frequency of 
the Hopf bifurcations in and e (data not shown). 




FIG. 2 Single patch behavior in the absence of cross immunity: (a) ADE value where the Hopf 
bifurcation occurs versus contact rate 0. (b) Period of the periodic orbit versus for = 1.903 . 
e = 0, v = 0, fi = 0.02, and other parameters as in Table HI 
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IV. COUPLED PATCHES 



We next examine the effect of coupling between distinct regions on the dynamics previ- 
ously observed in the single patch model. In particular, we consider how migration between 
nonidentical patches affects the stability of the steady state. We investigate the dynamics 
of the coupled systems by nu merically integ r ating Eqs. [T]{5] and by tracking the bifurcations 
using a continuation routine (IDoedel et all 119971 ). As mentioned in the previous Section, 
the model has a Hopf bifurcation at critical values of the bifurcation parameters <p, the ADE 
factor, and e, the cross immunity strength. We study the effect of coupling and asymmetry 
on the stability by observing their effect on the location of the Hopf bifurcations. We con- 
sider cross immunity and ADE separately; that is, we analyze the system either with ADE 
and no cross immunity or with cross immunity and no ADE. Including both cross immunity 
and ADE together leads to qualitatively similar behavior to that reported here, at least 
locally when the asymmetry is not too large. 

We proceed by fixing the patch-specific parameters in patch 2 and varying a social pa- 
rameter (Pi or fii) in patch 1 to observe the effect of increasing asymmetry on the dynamics. 
As previously mentioned, changing the social parameters modifies the natural frequency of 
the system. For each value of asymmetry, we look for the critical points at which the coupled 
system loses stability. We also consider several migration rates v. 

The effect of asymmetry in the contact rate is shown in Figs. [3]^a) and [3](b) for two 
migration rates, namely v = 0.02 (black) and v = 0.05 (gray). The value of the critical 
parameter e (Fig. EJ^a)) or </> (Fig. E](b)) at which the Hopf bifurcation occurs is plotted 
against the varying contact rate. The symmetric case is at the bottom of the curve. We see 
that two identical patches have the same bifurcation point as a single, well mixed population 
even when migration is present (cf. Fig. [2^). However, a striking difference in the dynamics 
appears when we make the two patches weakly asymmetric. The values of </> and e at which 
the Hopf bifurcation occurs are dramatically different from the symmetric case. The steady 
state stability persis ts well into the parameter regime where a single patch would display 



chaotic oscillations (Bi anco et al. 



2009 



Billings et al.l . 120071 ). The location of the single 



patch Hopf point does not depend strongly on /?, and decreasing the contact rate is actually 
destabilizing (Figure^), so the stabilization observed here is clearly a result of the coupling 
between asymmetric systems. Slightly increasing the migration rate (from v = 0.02 to 
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v = 0.05 in Figure [3]) further increases the stability of the system. 

Similar results are obtained if the asymmetry occurs in the other social parameter, the 
birth rate. Again, for asymmetric patches either stronger cross immunity or stronger ADE is 
needed to destabilize the steady state through a Hopf bifurcation (figure in Supplementary 
Material). 
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Unstable steady state 
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V 



Stable steady state 



(b) 



Unstable steady state 



TV 
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FIG. 3 Critical values of the parameters (a) e and (b) cj> at which the Hopf bifurcation occurs for 
coupled patches, as a function of the contact rate in patch 1, for two values of the migration rate 
v = 0.02 (solid black) and v = 0.05 (dashed gray). The contact rate in patch 2 is fixed at /?2 = 200. 
In (a), <p = 1 (no ADE). In (b), e = (no cross immunity). ^1=^-2 = 0.02 and other parameters 
are as in Tabled 



V. ANALYTICAL MOTIVATION 

We now motivate the results of the preceding Section using a simple lower dimensional 
model and show that the qualitative features depend only on the bifurcation structure and 
characteristic frequencies rather than on details of the epidemic model. In this Section we 
study the general behavior of two coupled systems, each displaying a Hopf bifurcation, but 
with different characteristic frequencies. 



The generic form for a Hopf bifurcation is the following, in polar coordinates (jStrogatz , 



2001) 

(6) 

9 =uj 

A Hopf bifurcation occurs at A = 0, where periodic oscillations with frequency u are ob- 
served. We now couple two systems of the sort in Eq. [61 each with a potentially different 
frequency u q , via linear migration terms with migration rate v. To lowest order (not dis- 
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playing cubic terms), the coupled system in cartesian coordinates is 



Xx\ — uJiyi — vx\ + vx 2 
Xyi + uixi - vy x + vy 2 

\X 2 - UJ 2 y 2 - VX 2 + VX X 

Xy 2 + uj 2 x 2 - vy 2 + vy x 



(7) 



Stability is determined by evaluating the Jacobian of Eq. [7] at the steady state 
(xi,yi,x 2 ,y 2 ) = 0. At the Hopf bifurcation, the real part of the largest eigenvalue crosses 
zero, with nonzero imaginary part. We obtain the Hopf point directly by setting the real 
part of the eigenvalue to zero and solving for the critical value of A at the bifurcation. Details 
may be found in the Supplementary Material. The Hopf bifurcation occurs at 



In Fig. H] we show the location of the Hopf bifurcation given by Eq. [8] as a function of 
the asymmetry between the two systems for v = 0.05. Comparison with Fig. [3] (and Fig. 1 
of the Supplementary Material) shows the qualitative agreement of the theoretical results 
for the lower dimensional system with the full multistrain system. Increasing the migration 
rate v increases the value of the bifurcation parameter to which the Hopf point saturates 
(not shown), as in Fig. [3j It is also worth noticing that, in the case of identical frequencies 
u)\ = u 2 , the Hopf bifurcation occurs at A c = 0, the location in the absence of coupling. This 
is consistent with the numerical results for the multistrain system in the case of symmetric 
patches. 

VI. CONCLUSIONS AND DISCUSSION 

We have studied the endemic steady state stability properties for a multistrain epidemic 
model on two migration-coupled patches. Interactions between strains in the model were 
governed by temporary partial cross immunity and antibody-dependent enhancement. In 
the absence of coupling, the system displayed Hopf bifurcations in two epidemic parameters. 
Coupling between patches with non-identical parameters, which gave them non-identical 
characteristic frequencies of oscillation, was shown to shift the Hopf bifurcations, stabilizing 
the steady state. This behavior was observed for the Hopf bifurcation obtained by sweeping 




(8) 
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Stabl9 steady state 

" '°3.8 3.85 3.9 3.95 4 4.05 4.1 4.15 4.2 



FIG. 4 Critical values of the parameter A at which the Hopf bifurcation occurs for the coupled 
generic Hopf bifurcations, as a function of the frequency Here the parameter U02 has been kept 
fixed at u>2 = 4. v = 0.05. 



the cross immunity in the absence of ADE and the bifurcation obtained by sweeping the 
ADE in the absence of cross immunity. It occurred for asymmetry in either of our two social 
parameters, the birth rate and the contact rate. 

To motivate this result, we diffusively coupled two low dimensional Hopf bifurcations 
with different characteristic frequencies and analyzed the stability of the steady state. We 
again saw that coupling between asymmetric systems led to stabilization. This indicates 
that the stabilization in the epidemic model is a result of the underlying dynamics, rather 
than the details of the model. The stabilization may occur because when the frequencies 
of oscillation in two coupled systems are different, the oscillations tend to cancel each other 
because of phase differences. This topic will be studied in more detail in a future work. 

Bifurcations from steady state to oscillatory behavior can be associated with an increased 
number of infection cases, particularly if chaoti c oscillations occur, as in previous dengue 



models (iBianco et al. 



2009 



Billings et all I2007I ). Our results suggest that if control strate- 



gies in one region are able to generate enough asymmetry, this could lead to a stabilization 
of the outbreaks, which would have a positive effect on adjacent regions. Asymmetry could 
be generated in the effect ive contact ra te by mosquito control, which could include reducing 



mosquito bre eding sites (jSlose 



development (iBarbazan et al. 



1986I ) . or through new genetic controls which are under 



2008I ). Asymmetry in the birth rate could be generated by 



lowering the effective birth rate through vaccination of new susceptibles once a vaccine is 
available. However, because the bifurcation point saturates rather than increasing indef- 
initely as asymmetry increases, such a strategy would be successful only if the epidemic 
parameters in the real system are moderately close to the bifurcation. (Extreme asymme- 
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try may even be destabilizing, so this strategy works best locally when the asymmetry is 
not too strong.) In addition, the role of seasonality in exciting oscillations should not be ig- 
nored. Seasonal variat ions in the contact rate have been included in previous dengue models 



(jSchwartz et al 



20051 ). The interplay of seasonality and coupling is a topic for future study. 



The results of the present work may also be useful in lower bound estimation of the 
parameter values for ADE and cross immunity. ADE especially is difficult to measure in 
vivo and must be estimated by other means. Since recent publ ications have su g geste d that 



epidemiological data from Thailand show chaotic outbreaks (jSchwartz et al. 



20051 ). and 



given the certainty of human migration and asymmetry between adjacent regions of the 
country, it is possible that the actual parameter values for ADE and cross immunity are 
higher than the ones estimated by studying a single well mixed population. 

Finally, the work discussed here shows a potential effect of human movement between 
heterogeneous regions. As spatial effects are further studied in epidemic models, it remains to 
be seen how this phenomenon will extend to more complicated spatial geometries, including 
more patches and perhaps non-symmetric coupling terms. Because the migration-induced 
stabilization depends only on the existence of a Hopf bifurcation in the model, it is expected 



that the stabiliza tion will be obse r ved in other population mode 
bifurcations (e.g., 



Fussmann et al. 



(I2000h : iGreenhakh et all (120041 ) I 



s that also contain Hopf 
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FIG. 5 Critical values of the parameters (a) e and (b) cj> at which the Hopf bifurcation occurs for 
coupled patches, as a function of the birth/death rate in patch 1, for v = 0.02. The birth/death 
rate in patch 2 is fixed at [12 = 0.02. In (a), = 1 (no ADE). In (b), e = (no cross immunity). 
(3\ = /?2 = 200 and other parameters are as in Table I of the paper. 



APPENDIX A: Asymmetric birth rates 

The effect of asymmetric birth rates is shown in Fig. of this supplementary material. 
The birth rate in patch 2 was held fixed while that in patch 1 was varied. For asymmetric 
patches either stronger cross immunity or stronger ADE is needed to destabilize the steady 
state through a Hopf bifurcation, compared to the symmetric patch case. 



APPENDIX B: Details of Hopf bifurcation analysis 

We evaluate the Jacobian for the ODE system given in Eq. 7 of the paper at the steady 
state {xi, yi,x 2 , 1/2} = 0. The roots of the characteristic polynomial f(t) of the Jacobian are 
the eigenvalues {t}. The four eigenvalues are 



ti/2 = \ - v + ]^\J Av 2 - 2{u\ + to 2 ) ± 2(cg>i + lo 2 ) \/iuh - to 2 ) 2 - 4^ 2 (Bl) 
t m = X-u-^Au 2 -2(ul + ui)±2(u 1 + u 2 ) v ^H~-uj2) 2 -4:U 2 , (B2) 

but it is not easy to see from these expressions where the Hopf bifurcation occurs. Instead, 
we solve directly for the Hopf point. 

At a Hopf bifurcation, a pair of eigenvalues has zero real part. Thus to obtain the Hopf 
point, we can set t = ifx, where /i is real. If ifi is an eigenvalue, we require f(ifi) = 0, so 

»[/(</*)] = (B3) 
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and 

9f[/(Vi)] = (B4) 

If A ^ v, Eq. [B4] has nonzero roots fx = ±|-y/4A 2 — 8Xu + 2uj\ + 2o>f. When this 
is substituted into Eq. IB31 we obtain four potential roots for A at the Hopf point. Two 
are complex, which is unphysical and can be ignored. The one corresponding to the Hopf 
bifurcation is A c = v — | a/4z/ 2 — (ui — U2) 2 ■ (The fourth root is larger and corresponds to 
loss of stability of the second pair of eigenvalues.) This A c is real and thus gives the Hopf 
point location when \u\ — t^l < 2za When uo\ — a>2 7^ and v > 0, A c is positive, indicating 
stabilization due to the asymmetry and coupling. 

When A = v, Eq. lB4l is always satisfied. In that case, we must turn to Eq. IB3l to determine 
/i. Eq. IB3I has four roots, fi = ±|(^i + 0J2) ± §\/ (^1 — — 4i> 2 . jj, must be real by our 
assumption, and this occurs when |a>i — uj 2 \ > 2u, which is precisely the case not covered by 
the above result. Thus when \u\ — lu 2 \ > 2u, the Hopf point is at A c = v. 
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